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Abstract 



We investigate the dynamics of the chiral transition in an expanding quark-anti- 
quark plasma. The calculations are made within a linear cr-model with explicit quark 
and antiquark degrees of freedom. We solve numerically the classical equations of 
motion for chiral fields coupled to the fluid dynamical equations for the plasma. Fast 
initial growth and strong oscillations of the chiral field and strong amplification of long 
wavelength modes of the pion field are observed in the course of the chiral transition. 

ON 

Introduction. - It is commonly believed that color deconfinement and chiral symmetry 
restoration take place at early stages of ultra-relativistic heavy-ion collisions. At interme- 
P-i! diate stages of the reaction the quark-gluon plasma may be formed and evolve through the 
states close to thermodynamical equilibrium. However, at later stages of the expansion 
the transition to the hadronic phase with broken chiral symmetry should take place. The 
breakdown of chiral symmetry will possibly lead to such interesting phenomena as formation 
of disoriented chiral condensates (DCCs) and classical pion fields, as well as clustering of 
quarks and antiquarks. These phenomena were studied recently in many publications [1-15], 
using QCD motivated effective models, such as the linear and nonlinear a-models and the 
Nambu-Jona-Lasinio (NJL) model. Of course, these models have some significant shortcom- 
ings, e.g. they do not possess color confinement and the NJL-model is non-renormalizable. 
The key point is, however, that these models obey the same chiral symmetry as the QCD 
Lagrangian. 

In most applications of the u-model the quark degrees of freedom are disregarded (see 
e.g. [2-9,15]). The inclusion of quarks |10, [IT], 0] makes it possible to study the hadronization 



process. In this letter, we construct a relativistic mean field fluid dynamical model based on 
the linear a-model assuming that the quarks and antiquarks are in local thermodynamical 
equilibrium before, during, and after the chiral transition. This is an idealization which 
certainly cannot be true at later stages of the expansion. On the other hand, in [|1(], [TTJ, [12] 
we have considered another limiting case assuming a collisionless approximation for quarks 
and antiquarks. Below, we assume that the expansion can be well approximated by a 
spherical scaling solution at this stage. 
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The linear a -model. - To describe the complex dynamics of the chiral transition we adapt 
the linear cr-model where quarks and antiquarks interact with the classical a and 7? fields. 
Together, the scalar field a and the pion field 7? = (tti, ^2,^3) form a chiral field $ = (cr, 7?). 
The equations of motion obtained from the linear cr-model Lagrangian are: 

d^a(x) + X 2 [o- 2 (x) + t? 2 (x) - v 2 ]a(x) -H= -gp,(x), 

d^n(x) + \ 2 [a 2 (x) + vf 2 (x) - v 2 ]n(x) = -gp p {x), (1) 

where p s = (qq) and p p = z (575777) are the scalar and pseudoscalar densities generated by 
valence u and d quarks. They are determined by the interaction of quarks with meson fields. 
As shown in |l(J these densities can be expressed as p s (x) = ga{x)a{x) and p p = gn(x)a(x), 



with the scaling function 



a ( x ) = / / = d P = f K^. P) + nqfa P)]- ( 2 ) 



Here ^ = 2 x 2 x 3 = 12 is the degeneracy factor of quarks and m{x) is the constituent 
quark mass which is determined self-consistently through the meson fields, 

m 2 (x) =g 2 \a 2 (x) + vf 2 (x)]. (3) 

In eq.(0) n q and n-q are the quark and antiquark occupation numbers, which in thermal 
equilibrium are given by the Fermi-Dirac distributions. This linear a-model is approximately 
invariant under chiral SUl{2) ® SUr(2) transformations if the explicit symmetry breaking 
term Ha is small. The parameters are chosen in such a way that the chiral symmetry is 
spontaneously broken in the vacuum so that the expectation values of the meson fields are 
(a) = f n and (7?) = 0, where f n = 93 MeV is the pion decay constant. The constant H is 
fixed by the PCAC relation that gives H = f n m 2 , where = 138 MeV is the pion mass. 
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Then one finds v 2 = f 2 — ^f. X 2 is determined by the sigma mass, m 2 = 2X 2 f 2 + m 2 r , which 
we set to 600 MeV yielding A 2 ~ 20. The coupling constant g is fixed by the requirement 
that the constituent quark mass in vacuum, m vac = gf n , is about 1/3 of the nucleon mass, 
giving g w 3.3. With these parameters a chiral phase transition is predicted at T c w 132 
MeV [|TTJ. 

Fluid dynamical equations. - The relativistic Vlasov equation consistent with the linear 
cr-model can be written as: 



d ldm 2 (x) d 



dxn 2 dxn dp. 



f(x,p)=I oM , (4) 



where I co u is the collision integral and f(x,p) the scalar part of the fermion distribution 
function. Applying the moment expansion method and assuming local thermal equilibrium 
(Icon = 0) one arrives at the fluid dynamical equations of the form: 



dpT"* + a -d,m 2 = 0, (5) 

The energy-momentum tensor can be expressed in terms of the rest frame energy density 8 
and pressure P of quarks as: 

T iw = (£ + p) u »u v - Pgv. (6) 
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where is the collective 4-velocity. £ and P read: 



and 



Here m is given by eq.(|3]), and n q (p) and n-q(p) are the quark and antiquark occupation 
numbers: 



n q {p) 



exp 



y/p 2 + m 2 - u\ 1 ' 

— f — + 



, Uqip) = n g {p,fi^ -fx). (9) 



In this way ([!]), ([3]) and ([5]) constitute a self-consistent set of differential equations which 
can be solved numerically for meson fields and quark distributions. Below we put the quark 
chemical potential to zero, /j,=0 (baryon-free plasma). 

In this letter we assume spherical scaling expansion: = x M fr and r = Vt 2 - W. All 
quantities now depend on r only: cr(r), 7t(t), m(r), £{m{r) 1 T{r)) and P(m(r),T(r)). The 
d'Alembertian appearing in eqs.(^) is reduced to <9 M <9 M — > + - 4-. Moreover, eq.(|5|) 
reads: 

dS £ + P a dm 2 , . 

By introducing the dimensionless temperature 6 = — this equation can be represented as 

d£ L _d£_dm 2 _ d_£_dB_ 
dr dm 2 dr 39 dr 

6 can now be found by solving an equation of the form 

IdO _ dlnm 2 F 2 {9) 

T 2 Tt- Fi{9) ^T + —' (12) 

where F\{0) and ^(0) are known functions. One can easily find analytical solutions in two 
limiting cases. In the massless (ultrarelativistic) limit (9 = ^ S> 1), equation (|l|) reduces 

to 

1 d9 1 dO 9 

which yields 

T = ^T . (14) 

T 

In the Boltzmann (nonrelativistic) limit (6 = — <C 1), on the other hand, one obtains a 
logarithmic behavior, 

1 cZ0 3 , 1 11 , r 

¥T T = -^-2 0) ^r --- 3h ^ < 15) 

which can be rearranged to give an expression for T as a function of r, 

nr) = J^T- (16) 

Initial conditions.- In the chirally symmetric phase, the mean values of meson fields 
are close to zero. We take initial values of fields and their derivatives from the Gaussian 
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T = 300 MeV 



T = 500 MeV 




Figure 1: Meson fields in units of f w as functions of proper time shown for two different initial 
temperatures. The initial conditions for the left figure are Tq = 3 fm/c, T(r ) = 300 MeV, 
$(t )=(15, 15,0,0) MeV and (<9$/<9r) T=TO = ( 15000, 15000,0,0) MeV/fm. For the right figure, 
t = 2 fm/c, T(r ) = 500 MeV, $(r ) = (10,10,0,0) MeV and (d$/dr) r=T0 = (20000,20000,0,0) 
MeV/fm. 
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distributions around the mean values. The probability distribution of field fluctuations 5a 
and 5tt can be expressed through the change of the thermodynamical potential Q: 



W(8a, Sir) ~ exp 



AQV 



(17) 



T 

To second order in field fluctuations, the expansion of the thermodynamical potential yields 




T 2 

2 1 \2 p2 

9 "qj2 ~ X f* 



2 = T2" (T ~ Tc) 2 • (18) 



Thus, the initial fluctuations of the chiral field can be chosen within a sphere in isospace 
with radius equal to the variance of the Gaussian distribution, 

@var — 5 77tv) t>9\t / " V / 

# 2 z/ 9 (T 2 - Tl)V 

Numerical results.- We ran the simulation for several sets of randomly oriented initial 
conditions taken at the initial time r and temperature T = T(r ). The time r c at which the 
temperature drops below the critical temperature T c depends on the initial values of r and 
T . We have observed the following common trends. At low values of r < r c , i.e. at high 
temperatures corresponding to the symmetric phase, the a and the 7Tj fields are essentially 
identical, and they therefore oscillate with the same frequency and more or less the same 
amplitude (depending on how the energy is distributed between da/dr and dir^/dr). 

However, when the temperature approaches its critical value around r c = 7 fm/c and 
chiral transition starts, the field dynamics change. The a field rises rapidly, within 2-3 
fm/c, and comes close to the vacuum value. After the transition, it oscillates around its 
vacuum value, The 7Tj components continue to oscillate around zero, but with decreasing 
frequency. Whereas the a field oscillates with approximately the same period (~ 2 fm/c) 
as before the transition, the 7Tj oscillations slow down and show longer periods (~ 9 fm/c). 
This is easy to understand because in normal vacuum the pion mass is about 4 times smaller 
than the sigma mass. With the given initial values, the pion fields die out rather quickly 
(particularly for high initial temperatures). 

In Figure § we show the temperature evolution for several initial temperatures, mainly 
to point out the uniformity of the evolution, which is to be expected when in each case 
we choose r so that T = T c = 132 MeV is reached at around r c = 7 fm/c. Initially, the 
temperature oscillates around a decreasing curve which falls off like in the massless limit, 



as in eq.(|Tj). But after the transition it slows down as the constituent mass rises, without 
quite reaching the behavior of the Boltzmann limit. The left figure depicts the temperature 
curves for the two initial conditions shown in figure [l|. 

The right part of Figure ^ shows the evolution from several different initial states. The 
curves in the right frame don't exhibit the numerical dips of the T =500 MeV curve on the 
left because we have distributed the initial field strength and energy more evenly among all 
four chiral degrees of freedom. The curves overlap because we have chosen tq so that in each 
case the temperature reaches its critical value around r c = 7 fm/c. This figure emphasizes 
the universality of the evolution. 

Figure |3] shows the scalar density p s = ga(T)a(r) and the scaling function a(r) (as 
defined in eq.(^), but now as a function of r), for the two scenarios considered in Figure 
[[]. For an initial temperature T = 500 MeV, a(r) starts out at about twice as large a 
value as for T = 300 MeV, but once it reaches the starting point of the T = 300 MeV 
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Figure 2: Temperature as function of proper time for different initial temperatures. Initial 
conditions of left figure as in Figure [I]. The sharp dips are purely numerical artifacts. On the 
right, initial field fluctuations are distributed more evenly between the four chiral degrees 
of freedom. 
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Figure 3: Scalar density p s (r) and scaling function a(r). Initial conditions as in figure |I[ 
The sharp dips in the scaling function are purely numerical artifacts. 
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Figure 4: The power of the modes: k = 10,80,150 MeV translated into the (t,r) space. 
They follow the rule of the smaller the mode the larger the amplitude 
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curve, it has decreased and they evolve close to each other. From Figure || one can conclude 
that the scalar density of quarks and antiquarks drops out very rapidly after the transition 
temperarure. At later times the meson field dynamics is very little disturbed by the presence 
of quarks and antiquarks. 

We have calculated the power spectrum of the first component of the pion field for the 
modes with wave numbers k = 10, 80, 150. In the present context the power spectrum was 
first discussed in [HI. It is defined as 



S{t,k) 



47T f R 

— — / rdr sm(kr)ir(t, r) 
k Jo 



(20) 



where R = \/t 2 — r 2 is the spatial radius of the system at time t. The initial conditions are 
the same as in Figure [I]. In the case of an initial temperature of 300 MeV, one clearly observes 
the strong amplification of the low momentum mode as found in many other works |4], |5], [J . 
But we see that in the present case, partly due to the 3D expansion, the amplification of 
the pion field is weaker than found previously. Also, the response of the system depends 
sensitively on the initial conditions. This is clearly demonstrated by the power spectrum for 
T = 500 MeV, where the response of the low momentum mode is reduced by an order of 
magnitude compared to the case of 300 MeV. Clearly the field fluctuations will grow as the 
system approaches the critical temperature, and then they will be further amplified in the 
course of the chiral transition. In the present schematic model this effect is not taken into 
account because we define the initial fluctuations deep inside the symmetric phase and then 
propagate them through the chiral transition. 

Summary.- We have investigated the dynamics of the chiral transition in expanding 
quark- ant iquark plasma. The dynamics of chiral symmetry breaking and the time evolution 
of the chiral field were studied within a linear cr-model in a relativistic mean field fluid 
dynamical approach. This approach made it possible to study the evolution of chiral fields 
even before the chiral transition. Fast initial growth and strong oscillations of the chiral 
field and strong amplification of long wavelength modes of the pion field were observed. We 
are working on 3D simulations and taking dissipation into account. 
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